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We construct a quantum random walk algorithm, based on the Dirac operator instead of the Lapla- 
cian. The algorithm explores multiple evolutionary branches by superposition of states, and does 
not require the coin toss instruction of classical randomised algorithms. We use this algorithm to 
search for a marked vertex on a hypercubic lattice in arbitrary dimensions. Our numerical and 
analytical results match the scaling behaviour of earlier algorithms that use a coin toss instruction. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

Random walks are a fundamental ingredient of non- 
deterministic algorithms [3| , and are used to tackle a wide 
variety of problems — from graph structures to Monte 
Carlo samplings. Such algorithms have many evolution- 
ary branches, which are explored probabilistically, to esti- 
mate the correct result. A classical computer can explore 
only one branch at a time, so the algorithm is executed 
several times, and the estimate of the final result is ex- 
tracted from the ensemble of individual executions by 
methods of probability theory. Such algorithms are typi- 
cally represented using graphs, with vertices denoting the 
states and the edges denoting the evolutionary routes. A 
particular evolution corresponds to a specific walk on the 
graph, and the final result is obtained by combining the 
results for many different walks. To ensure that differ- 
ent evolutionary branches are explored in different execu- 
tions, one needs non-deterministic instructions, and they 
are provided in the form of random numbers. A coin toss 
is the simplest example of a random number generator, 
and it is included in the instruction set for a probabilistic 
Turing machine. 

A quantum computer can explore multiple branches 
of a non-deterministic algorithm in a single attempt, by 
using a clever superposition of states. The probabilistic 
result can then be arrived at by interference of ampli- 
tudes corresponding to different branches. Thus as long 
as the means to construct a variety of superposed states 
exist, there is no a priori reason to include a coin toss 
as an instruction for a (probabilistic) quantum Turing 
machine. 

In what follows, we construct a quantum random walk 
on a hypercubic lattice in arbitrary dimensions without 
using a coin toss instruction, analyse its properties, and 
use it to find a marked vertex on the lattice. More details 
are available in Refs.0,Q. 
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II. DIFFUSION 

A random walk is a diffusion process, commonly de- 
scribed using the Laplacian operator in the continuum. 
To construct a discrete quantum random walk, we must 
discretise the diffusion process using evolution operators 
that are both unitary and ultra-local (an ultra-local op- 
erator vanishes outside a finite range). 

On a periodic lattice, the spatial modes are charac- 
terised by discrete wave vectors k. Quantum diffusion 
then depends on the energy of these modes according 
to U{k,t) ~ cxp{—iE{k)t). The lowest energy mode, 
k = 0, corresponding to a uniform distribution, is an 
eigenstate of the diffusion operator and does not propa- 
gate. The slowest propagating modes are the ones with 
smallest nonzero |fc|. The classical Laplacian operator 

gives E{k) (X |A?|2 g, which translates to the character- 
istic Brownian motion signature, spread (n)ims « \/t. 
There is an alternative in quantum theory — instead of 
the non-relativistic Schodinger equation based on the 
Laplacian operator V^, one can use the relativistic Dirac 
equation based on the operator y. The Dirac operator 
gives E{k) cx |fc|, with the associated signature, spread 
(n)rins t. Clearly the Dirac operator, with its faster dif- 
fusion of the slowest modes compared to the Laplacian, 
is the operator of choice for constructing faster diffusion 
based quantum algorithms. 

An automatic consequence of the Dirac operator is 
the appearance of an additional internal degree of free- 
dom corresponding to spin, whereby the quantum state 
is described by a multi-component spinor. These spinor 
components were identified with the states of a coin in 
Refs.0,0], with the coin evolution rule guiding the quan- 
tum diffusion process. While this is the correct proce- 
dure in the continuum theory, another option is available 
for a lattice theory, i.e. staggered fermions 0- In this 
approach, the spinor degrees of freedom are spread out 
over an elementary hypercube, location dependent signs 
appear in the evolution operator, and translational in- 
variance exists in steps of 2 instead of 1. We follow this 
approach to construct, a quantum diffusion process on a 
hypercubic lattice, without a coin toss instruction. 
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The free particle Dirac Hamiltonian in d-space dimen- 
sions is 



W = UeUo 



^frcc — —'^^ • V + f3m 



(1) 



On a hypercubic lattice, the simplest discretisation of the 
derivative operator is 



1 



ykf{x)^-[fix + k)-f{x~k)] 



(2) 



Then the anticommuting matrices a, (3 can be spin- 
diagonalised to the location dependent signs 



fe-i 



/3 



(3) 



Even when the Hamiltonian H is ultra-local (i.e. has 
a finite range), the evolution operator U = exp{—iHt) 
is not. To make the evolution operator ultra-local, we 
break up H in to block-diagonal Hermitian parts, and 
then exponentiate each part separately. Partitioning of 
H in to two parts (which we label "odd" and "even") is 
sufficient for this purpose . This partition is illustrated 
in Fig.f for d = 1 and d = 2. Each part contains all the 
vertices but only half of the links attached to each vertex. 
Consequently, each link appears in only one of the two 
parts, and can be associated with a term in H providing 
propagation along it, i.e. 



free 



— Ho + He 



(4) 



The Hamiltonian is thus divided in to a set of non- 
overlapping blocks that can be exponentiated exactly. 
Each block is an elementary hypercube on the lattice, 
and the block matrices are of size 2^ x 2"^ in d dimensions. 
The ultra-local quantum random walk on the lattice then 
evolves the amplitude distribution according to 



V'(f;0 = W*ij{x; 0) 



(5) 



(a) 



(b) 






FIG. 1: Partitioning of lattice in to odd and even parts for 
ultra-local discrete quantum walk: (a) for d = 1, (b) for d = 2. 



(6) 



Each block of the unitary matrices ?7o(e) mixes the ampli- 
tudes of vertices belonging to a single elementary hyper- 
cube, and the amplitude distribution spreads because the 
two alternating matrices do not commute. The random 
walk operator W is translationally invariant in steps of 
2, along each coordinate direction. 



III. QUANTUM RANDOM WALK ON A LINE 
A. Construction 

To explicitly illustrate the above described procedure, 
let us consider the random walk on a line, with the al- 
lowed positions labeled by integers. 

The simplest translation invariant ultra-local discreti- 
sation of the Laplacian operator is 



H\n) cx [- |n- 1) +2|n) - |n-|- 1)] 



(7) 



One may search for ultra-local translationally invariant 
unitary evolution operators using the ansatz 



U\n) = a\n - 1) 



b\n) + c\n + 1) 



(8) 



but then the orthogonality constraints between different 
rows of the unitary matrix make two of {a, 6, c} vanish, 
and one obtains a directed walk instead of a random walk. 
This problem can be bypassed, and an ultra-local unitary 
random walk can be constructed, by enlarging the Hilbcrt 
space with a quantum coin, e.g. 



U 



E[it)(Ti 



®|n+f)(n| + |i)(i|^|n-l)(n| . (9) 



This route brings its own set of caveats, due to quan- 
tum entanglement between the coin and the position de- 
grees of freedom. 

We follow an alternate route familiar to lattice field 
theorists JJ. ft has also been used to simulate quan- 
tum scattering with ultra-local operators |lO|. and to 
construct quantum cellular automata The starting 
point is the decomposition of the Laplacian operator in 
to its even and odd parts, H = Hp + Ho, 



H cx 



He OC 




-1 
f 




(10) 



(11) 
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H„ cx 
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-1 






-1 

1 



(12) 
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While H has the structure of a second derivative, its two 
parts, He and Ho, have the structure of a first derivative. 
The above decomposition is indeed reminiscent of the 
"square-root" one takes to go from the Laplacian to the 
Dirac operator. 

The two parts, H,, and Ho, are individually Hcrmitian. 
They are block-diagonal with a constant 2x2 matrix, and 
so they can be exponentiated while maintaining ultra- 
locality. The total evolution operator can therefore be 
easily truncated, without giving up either unitarity or 
ultra-locality. 



U{At) 



i{H^+H„)At 



iH„At iH^At 

e ' e ° 



0{iAtf)il3) 
U,{At)Uo{M) + 0{{Mf) . 



The quantum random walk can now be generated using 
UelJo as the evolution operator for the amplitude distri- 
bution ^{n, t), 



(14) 



The fact that Ue and Uo do not commute with each other 
is enough for the quantum random walk to explore all 
possible states. The price paid for the above manip- 
ulation is that the evolution operator is translationally 
invariant along the line in steps of 2, instead of 1. 

The 2x2 matrix appearing in He and Ho is propor- 
tional to (1 — (Ti), and so its exponential will be of the 
form (cl -I- isCTi), |cp-t-|sp = l. A random walk should 
have at least two non-zero entries in each row of the evo- 
lution operator. Even though our random walk treats 
even and odd sites differently by construction, we can 
obtain an unbiased random walk, by choosing the 2x2 
blocks of Ue and Uo as -^Ql)- Furthermore, it is com- 
putationally more convenient to choose a basis where the 
unitary operators are all real. Performing a global phase 
transformation, |n) — > i"\n the 2x2 blocks of Ue 

and Uo become -^(1 ± ia2)- The discrete quantum ran- 
dom walk then evolves the amplitude distribution accord- 
ing to 



UeUo\n) 



Uo\n) 
Ue\n) 

1 r 



1 

1 r 

V2 



|n)-(-l)"|n+(-l)") 
|n> + (-l)"|n-(-l)") 



(15) 
(16) 



|n-l) + |n)-|n+l) + |n + 2(-l)") . (17) 



It is instructive to realise that, with the above choice, 
the unbiased quantum random walk represents the path 
integral for a relativistic particle with \p\ = m. Its speed 



(in units of speed of hght) is then |w| = 1/a/2. The 
directed walk, with the 2x2 block matrix [/ cx (T2, cor- 
responds to \v\ = 1, and the stationary limit U = 1 
corresponds to w = 0. 



B. Analysis 

It is straightforward to analyse the properties of the 
walk in Ea. (|17|l using the Fourier transform: 



dk 



>(fc,i) 



(18) 



(19) 



The evolution of the amplitude distribution in Fourier 
space is easily obtained by splitting it in to its even and 
odd parts: 



— ( i'e 



, V(fc,i) = [M(fc)]*V(fc,0) 



M{k) 



cosfc —ismk 
-isinfc e~*'^ cosfc 



(20) 



(21) 



The unitary matrix M has the eigenvalues, A± = e^*'^'= 
(this ± sign label continues in all the results that follow) , 



A± = cos^ fc ± i sin k \/l + cos^ k , ujk = cos ^ (cos^ k) , 

(22) 

with the (unnormalised) eigenvectors. 



e± oc 



cos fc ^ V 1 + cos^ fc 
1 



cos fc =F Vl -I- cos^ fc 



(23) 



The evolution of amplitude distribution then follows 

V^(fc,t) = e™'=*V^+(fc,0) + e-™''*V^_(fc,0) , (24) 

where i^±{k,Qi) are the projections of the initial ampli- 
tude distribution along e±. The amplitude distribution 
in the position space is given by the inverse Fourier trans- 
form of tp{k,t). While we are unable to evaluate it ex- 
actly, many properties of the quantum random walk can 
be extracted numerically as well as by suitable approxi- 
mations. 

A walk starting at the origin satisfies -00(^^,0) = 6n,o- 
This walk is asymmetric because our definitions treat 
even and odd sites differently. We can construct a 
symmetric walk, using the initial condition '0s(n, 0) = 
{Sn.o + iSn,i)/V^- The resultant probability distribution 
is then symmetric under n (1 — n). (Real and imag- 
inary components of the amplitude distribution evolve 
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independently because we have chosen the evolution op- 
erator to be real.) For both these initial conditions, by 
construction, the quantum random walk remains within 
the interval [—2t + 1, 2t] after t time steps. 

The escape probability of the quantum random walk 
can be calculated by introducing a fully absorbing wall, 
say between n = and n = — 1. Mathematically, this 
absorbing wall amounts to a projection operator for n > 
0. The unabsorbed part of the walk is given by 



1 



(25) 



UeUo Hn,t) - -5„,_i(V'(0,t) + V(l,t)) 



with the absorption probability, 

Pabs(t) = l-^|V(",i)l' 



(26) 



n>0 



All these variations in initial and boundary conditions 
are easy to implement numerically, and examples are 
shown in Fig. 2. We have used such simulations to study 
various properties of the quantum random walk. 

For large t, a good approximation to the probability 
distributions can be obtained by the stationary phase 
method 0, l^. The smoothed probability distribution 
for the symmetric walk, obtained by replacing the highly 
oscillatory terms by their mean values, is 



s I smooth 



(27) 



(Here, the n <-> (1 — n) symmetry can be restored by 
replacing n by (n— i).) As shown in the top part of Fig. 2, 
it represents the average behavior of the distribution very 
well. Its low order moments are easily calculated to be. 



\/2t 



s I smooth 



dn 



\/2t 



-V2t 
V2t 



\n\ ■ lip: 



's I smooth 



.2,.,, ,2 , 



s I smooth 



1 , (28) 
t , (29) 
2(2 - V2)t^ . (30) 



Results 



The following properties of the quantum random walk 
are easily deduced 

• The probability distribution is double-peaked with 
maxima approximately at ±\/2t. The distribution falls 
off steeply beyond the peaks, while it is rather flat in 
the region between the peaks. With increasing t, the 
peaks become more pronounced, because the height of 
the peaks decreases more slowly than that for the flat 
region. The location of the peaks is in accordance with 



the propagation speed, \v\ = l/-\/2, once we take in to 
account the fact that a single step of our walk is a prod- 
uct of two nearest neighbor operators, Ue and Uo- 

• The size of the tail of the amplitude distribution is 
limited by (ei)^^ ~ t~^^^, which gives An> = A(et) — 
0{t^^^). On the inner side, the width of the peaks is 
governed by Iw^'tl-^/^ - f-i/^. For \n\ = {V^ - S)t, this 
gives An< = A{St) — 0{t^/'^). The peaks therefore make 
a negligible contribution to the probability distribution, 

0(^-1/3). 

• Rapid oscillations contribute to the probability distri- 
bution (and hence to its moments) only at subleading 
order. They can be safely ignored in an asymptotic anal- 
ysis, retaining only the smooth part of the probability 
distribution. 

• The quantum random walk spreads linearly in time, 
with a speed smaller by a factor of \/2 compared to a 
directed walk. This speed is a measure of its mixing 
behavior and hitting probability. The probability distri- 
bution is qualitatively similar to a uniform distribution 
over the interval [—^/2t, \/2t]. In particular, the m'^ mo- 
ment of the probability distribution is proportional to 
t™. This behaviour is in sharp contrast to that of the 
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FIG. 2: Probability distribution after 32 time steps for the 
symmetric quantum random walk ifja- In the top figure, the 
dark curve denotes the smoothed distribution of Ea. H27|l . The 
bottom figure shows the distribution in the presence of an 
absorbing wall on the left side of n = 0. 
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classical random walk. The classical random walk pro- 
duces a binomial probability distribution, which in the 
symmetric case has a single peak centered at the origin 
and variance proportional to t. The linear spread in time 
of our quantum random walk is achieved even when tp 
has 50% probability to stay put at the same location at 
every step, as can be seen from Eas. (|15ll6ll . This means 
that our walk is more directed and less of a zigzag. 

• Above properties agree with those obtained in Refs.j^ 
IT^ for a quantum random walk with a coin-toss instruc- 
tion (extra factors of 2 appear in our results because of 
difference in our conventions), demonstrating that the 
coin offers no advantage in this particular set up. Essen- 
tially, we have absorbed the two states of the coin in to 
the even/odd site label at no extra cost. By making the 
coin states part of the position space, we have eliminated 
quantum entanglement between the coin and the position 
degrees of freedom completely — only superposition rep- 
resenting the amplitude distribution survives 14]. Such 
a reorganisation would be a tremendous advantage in any 
practical implementation of the quantum random walk, 
because quantum entanglement is highly fragile against 
environmental disturbances while mere superposition is 
much more stable. The cost for gaining this advantage 
is the loss of short distance homogeneity — translational 
invariance holds in steps of 2 instead of 1. 

• Comparison of the numerically evaluated probability 
distributions in Fig. 2, without and with the absorbing 
wall, shows that the absorbing wall disturbs the evolu- 
tion of the walk only marginally. The probability distri- 
bution in the region close to n = is depleted as antici- 
pated, while it is a bit of a surprise that the peak height 
near n — increases slightly. As a result, the escape 
speed from the wall is little higher than the spreading 
speed without the wall. Overall, the part of the quan- 
tum random walk going away from the absorbing wall 
just takes off at a constant speed, hardly ever returning 
to the starting point. Again, this behavior is in a sharp 
contrast to that of the classical random walk, which al- 
ways returns to the starting point, sooner or later. We 
also find that the first two time steps dominate absorp- 
tion, Ps.abs(i = 1) = 0.25 and Fs,abs(i = 2) = 0.375, with 
very little absorption later on. Asymptotically, the net 
absorption probability approaches Ps,abs(oo) « 0.4098 for 
the symmetric walk. This value is smaller than the corre- 
sponding result Pabs(oo) = 2/7r for the symmetric quan- 
tum random walk with a coin-toss instruction [T^. 



IV. QUANTUM RANDOM WALK 
ON A HYPERCUBIC LATTICE 



2-dim Lattice 



He\x,y) 
H\x,y) 



+ (-l)-+''|x,y+(-l)^)J , (31) 
{-lY\x-{-ir,y) 

+ (-l)-+^|x,y-(-l)«)] , (32) 

{Ho + H,)\x,y) 
1 r 

+ {-ir{\x,y + l)-\x,y-l))\ . (33) 



More explicitly, the 4x4 blocks of the Hamiltonian are: 



-1 -1 




-(/(g) a2 + a-2® (T3) 



(34) 
(35) 



where the column on the right denotes the vertices of 
the elementary square on which operates. Similarly, 
— —H^, when operating on the square with vertices 
{00,-10,0-1,-1-1}. Noting that = Hi = i/, the block- 
diagonal matrices are easily exponentiated to 



U, 



o(e) 



= cl 



= 1 



(36) 



The parameter c (or s) is to be tuned to achieve the 
fastest diffusion across the lattice. 

The quantum random walk with the Dirac operator 
spreads on a two-dimensional grid as illustrated in Fig. 3. 
The continuum Dirac Hamiltonian has exact rotational 
symmetry, and that survives to an extent even after dis- 
cretisation on a hypercubic lattice. After a point start, 
the random walk spreads essentially isotropically at dis- 
tances much larger than the lattice spacing, while the 
hypercubic symmetry governs the random walk pattern 
at shorter distances. Of course, the hypercubic symme- 
try would be exact for a d-dim random walk constructed 
as a tensor product of d one-dimensional random walks. 



.005 
. 004 
. 003 
. 002 
. 001 





/1 00 



Next let us consider the situation for d — 2. The par- 
titioned free Hamiltonian is given by 



Ho\x,y) 



i-ir\x+i-ir,y) 



FIG. 3: Probability distribution for the quantum random 
walk on a two-dimensional 128 x 128 grid after 61 time steps. 
A symmetric start, i.e. ip{x,y,Q) = {S^oSyo + iSxiSyi) / \/2, 
was used. Darker shades on the grey-scale indicate larger 
probability. 
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B. 3-dim Lattice 

Next let us look at the situation for d — 3. The parti- 
tioned free Hamiltonian is given by 



and = —H^ when operating on the hypercube with 
coordinates flipped in sign. The block-diagonal matrices 
satisfy = ^ J-^^ ^^'^ exponentiate to 



Ho\x,y,z) 



U, 



o(e) 



Cl 



IS 



H\x,y,z) 



{~ir\x + {-lf,y,z) 
-ir+y\x,y+{-l)y,z) (37) 

-ir+y+^\x,y,z + {-iy) 

(-inx-(-ir,y,z) 
-ir+y\x,y^{^l)y,z) (38) 

-ir+y+^\x,y,z-{-iy) , 

{Ho + H^)\x,y,z) 

\x + l,y,z) -\x- l,y,z) 
iri\x,y + l,z)-\x.,y-l,z) (39) 
+ {-ir+y{\x,y,z + l)-\x,y,z-l)) . 

More explicitly, the 8x8 blocks of the Hamiltonian are: 



Vd 



H, 



o(e) 



(43) 
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, (40) 



with the column on the right indicating the vertices 
of the elementary cube on which operates. Like- 
wise, — —H^, when operating on the elementary 
cube with vertices {000,-160,0-10,-1-10,00-1,-10-1,0-1-1, 
-1-1-1}. With Hl^ Hi ^ |/, the block-diagonal matri- 
ces exponentiate to 



Uo(e) = cl - is — Ho(^e) 



Icp-f ksP = l 



(41) 



V. SEARCH ON A HYPERCUBIC LATTICE 
USING THE DIRAC OPERATOR 

A. Strategy 

A clear advantage of quantum random walks is their 
linear spread in time, compared to square-root spread in 
time for classical random walks. So they are expected to 
be useful in problems requiring fast hitting times. Sev- 
eral examples of this nature have been explored in graph 
theoretical and sampling problems (see Refs.flsL llGj for 
reviews). Here we consider the particular case of using 
the quantum random walk to find a marked vertex on a 
hypercubic lattice (see also Refs.^^]). 

Consider a d-dim hypercubic lattice with N — L"^ 
vertices, one of which is marked. The quantum algo- 
rithmic strategy for the search process is to construct 
a Hamiltonian evolution, where the kinetic part of the 
Hamiltonian diffuses the amplitude distribution all over 
the lattice while the potential part of the Hamiltonian 
attracts the amplitude distribution towards the marked 
vertex 17|. The optimisation criterion is to concentrate 
the amplitude distribution towards the marked vertex as 
quickly as possible. In his algorithm, Grover constructed 
a global operator that allows diffusion from any vertex to 
any other vertex in just one step. Under different circum- 
stances, when diffusion is restricted to be ultra-local (i.e. 
one can only go from a vertex to its neighbours in one 
step), one must find an appropriate diffusion operator 
that provides fast propagation of spatial modes. Obvi- 
ously, the Dirac operator is better suited to this task than 
the Laplacian operator. 

To search for a marked vertex, say the origin, we need 
to attract the quantum random walk towards it. This 
can be accomplished by adding a potential to the free 
Hamiltonian, 



Again c (or s) is a parameter to be tuned to achieve the 
fastest diffusion across the lattice. 



C. d-dim Lattice 

We can now observe a pattern in the explicit results 
for d = 1, 2, 3 above. The 2'^ x 2"^ blocks of the Hamilto- 
nian can be written as sums of tensor products of Pauli 
matrices. As suggested by Eas. H35l4UII . 



(44) 



Exponentiation of this potential produces a phase change 
for the amplitude at the marked vertex. It is optimal to 
choose the magnitude of the potential to make the phase 
maximally different from 1, i.e. e~*^°'^ = —1, whereby 
the phase becomes a reflection operator (binary oracle). 



R = / - 2|0)(0| 



(45) 



(42) 



The search algorithm alternates between the diffusion 
and the reflection operators, yielding the evolution 



^{x-M,t2) = [T/F*ii?]*^V(^;0,0) . 



(46) 
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Here <2 is the number of oracle calls, and ti is the number 
of random walk steps between the oracle calls. Both have 
to be optimised, in addition to c and depending on the 
size and dimensionality of the lattice, to find the quickest 
solution to the search problem. 

Fastest search amounts to finding the shortest unitary 
evolution path between the initial state, typically chosen 
as the uniform superposition state |s) = '^xlx) /VN, 
and the marked state |0). This path is a circular arc 
(geodesic) from \s) to |0). With the random walk diffu- 
sion operator W, evolution of the state j?/") does not re- 
main restricted to the two-dimensional subspace formed 
|s) and |0). Thus to optimise our algorithm, we need to 
tune the parameters so as to 

(a) maximise the projection of the state lip) on to the 
two-dimensional \s) — |0) subspace, and 

(b) maximise the angle of rotation by the operator VF*^ R, 
for the projected component of |V') in the |s)-|0) subspace. 
We have explored this optimisation numerically. 




FIG. 4: Probability distribution for the quantum random 
walk search problem on a two-dimensional 64 x 64 grid, at 
the instance when the probability at the marked vertex at- 
tains its largest value. The number of random walk steps are 
ti = 1 (top) and ti = S (bottom). 



B. Numerical Results 

We carried out computer simulations of the quantum 
random walk search problem with a single marked ver- 
tex, for ci = 2, 3. The algorithm was optimised by tuning 
the parameters c and ti, so as to minimise the number 
of oracle calls ^2 required to find the marked vertex. The 
following is a summary of our observations: 

• An unbiased search starts with a uniform probability 
distribution over the whole lattice. Thereafter, the prob- 
ability at the marked vertex goes through periodic cycles 
of rise and fall as a function of time step. It is crucial 
to stop the algorithm at the right instance to find the 
marked vertex with a significant probability. 

• Our best results are obtained with c = 1/V2 and ti = 3. 
In this case, the probability at the marked vertex reaches 
its largest value, and t2 achieves its smallest value. With 
these parameters, the probability at the marked vertex 
shows a periodic sinusoidal behaviour, which persists for 
more than 30 cycles without any visible deviation. Also, 
apart from the uniform background, the probability dis- 
tribution shows a sharp single-point delta function at the 
marked vertex. These features indicate that the walk 
evolves largely in the two-dimensional subspace formed 
by the uniform state and the marked state. 

• For c < I/V2, the walk diffuses more slowly, and t2 
increases. For c > l/v^, the probability at the marked 
vertex loses its periodic sinusoidal behaviour, suggesting 
that the walk no longer remains confined to the two- 
dimensional subspace. For the optimal choice c = l/y/2, 
the probability of the walk remaining at the same vertex 
equals that for moving to a neighbouring vertex, which 
corresponds to the most efficient mixing between odd and 
even sublattices. 

• For ti < 3, the probability distribution spreads out in- 
stead of being a delta function at the marked vertex, as 
illustrated in Fig. 3. This decreases the peak probability 
at the marked vertex. Moreover, t2 increases. For ti > 3, 
the probability at the marked vertex loses its sinusoidal 
behaviour, again with a decrease in the peak probability. 
In both cases, the changes indicate that the walk is drift- 
ing out of the two-dimensional subspace. An appropriate 
choice of ti is thus crucial to keep the walk close to the 
two-dimensional subspace. 

• For the 2-dim walk, the largest probability at the 
marked vertex is predicted to be 0{1/ \ogN), which oc- 
curs after 0{y/N log N) time steps £,]. To make the 
marked vertex probability 0(1), an amplitude amplifica- 
tion procedure is required | l8| , and the overall search al- 
gorithm scales as 0{\/NlogN). Our numerical results, 
shown in Fig. 5, are consistent with these expectations. 
Simple fits provide the parametrisations: 

0(1/ log iV) ^ 2.12/log2Af , (47) 
0{^/N log N) — > 0.137 VA^logaiV . 

• For the walk in more than two dimensions, the largest 
probability at the marked vertex is predicted to be 0(1), 
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FIG. 5: The peak probability at the marked vertex for the quantum random walk search problem (left), and the number of 
time steps required to reach it (right), as a function of the database size. The dashed and continuous curves correspond to 
d = 2 and d = 3 respectively, c = 1/V2 and ti = 3 were used, with N ranging from 2^ to 2^*. 



which occurs after 0{y/N) time steps Our numer- 

ical results for the 3-dim walk, also displayed in Fig. 5, 
agree with these scaling rules. Simple fits provide the 
parametrisations: 

0(1) — > 0.0969 , 0{y/N) — > OMSVn . (48) 
These results demonstrate that our quantum random 



walk algorithm achieves the optimal scaling behaviour for 
the problem of finding a marked vertex on a hypercubic 
lattice. Thus our quantum random walk, based on the 
Dirac operator and not containing a coin toss instruc- 
tion, is no less effective in its diffusion properties than 
the earlier quantum random walks that use a coin toss 
instruction. 
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